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Context. One of the most dramatic events in the life of a low-mass star is the He flash, which takes place at the tip of the red giant 
branch and is followed by a series of secondary flashes before the star settles on the zero-age horizontal branch (ZAHB). Yet, no stars 
have ever been positively identified in this key phase in the life of a low-mass star (hereafter the "pre-ZAHB" phase). 
Aims. In this paper, we investigate the possibility that at least some pre-ZAHB stars may cross the instability strip, thus becoming 
variable stars whose properties might eventually lead to positive identifications. In particular, it has been suggested that some of the 
RR Lyrae stars with high period change rates (P) may in fact be pre-ZAHB stars. Here we present the first theoretical effort devoted 
to interpreting at least some of the high-P stars as pre-ZAHB pulsators. 

Methods. We constructed an extensive grid of evolutionary tracks using the Garching Stellar Evolution Code (GARSTEC) for a 
j_( . chemical composition appropriate to the case of the globular cluster M3 (NGC 5272), where a number of stars with high P values are 

found. We follow each star's pre-ZAHB evolution in detail, and compute the periods and period change rates for the stars lying inside 
' the instability strip, also producing pre-ZAHB Monte Carlo simulations that are appropriate for the case of M3. 

Results. Our results indicate that one should expect of order 1 pre-ZAHB star for every 60 or so bona-fide HB stars in M3. Among 
the pre-ZAHB stars, approximately 22% are expected to fall within the boundaries of the instability strip, presenting RR Lyrae-like 
pulsations. On average, these pre-ZAHB pulsators are expected to have longer periods than the bona-fide HB pulsators, and 76% 
of them are predicted to show negative P values. While the most likely P value for the pre-ZAHB variables is as -0.3 d/Myr, more 
t , extreme P values are also possible: 38% of the variables are predicted to have P < -0.8 d/Myr. 

Conclusions. It appears likely, therefore, that some - but certainly not all - of the RR Lyrae stars in M3 with high (absolute) P values 
■ are in fact pre-ZAHB pulsators in disguise. 

m ■ 

Key words, stars: evolution — stars: Hertzsprung-Russell diagram — stars: horizontal-branch — stars: oscillations — stars: variables: 
other — Galaxy: globular clusters: individual (M3 = NGC 5272) 
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C$ • ate helium cores as they continue their evolution accross the directly testing current models of pre-ZAHB evolution 

Hertzsprung-Russell diagram (HRD). As the mass of the core in- The reasons why it has proven so challenging to pinpoint 

creases due to the continued supply of He-rich material from the pre-ZAHB stars are twofold: first, the pre-ZAHB phase is very 

H-burning shell that surrounds the core, it also contracts and be- shorti compared to other major evolutionary phases in the life 

comes increasingly denser and hotter. Eventually, the conditions of a i ow . mass sta r. Second, pre-ZAHB stars are largely ex- 

in the core become ripe for the ignition of the He-burning triple- pec ted to overlap the loci occupied by asymptotic giant branch 

a reactions. Since the core at this point is electron-degenerate, (AGB), HB, and red giant branch (RGB) stars in observed color- 

the process does not take place quiescently, but rather in the form magnitude diagrams (CMD). 

of a thermonuclear runaway - the so-called He flash. The pri- Bearing these difficulties in mind , |ciiehn1 d2005l 12007b sug- 

mary He flash, which does not take place at the center proper gested an alternative approach for ±e detection of pre-ZAHB 

(due to the very efficient neutrino cooling there), is followed by namelyi through stdlar puLsations . In deed, some pre . 

a series of secondary flashes, each occurring increasingly closer ZAHB starfj are expected to rapidly cross the C epheid/RR Lyrae 

to the star's center, until degeneracy is lifted throughout the core instability strip on their route from the RGB tip t0 the ZAHB, 

and quiescent He burning finally commences - and a so-called thus becoming pu i sating stars a i ong the way. At least some of 

horizontal branch (HB) star is born (see | Catelan| | 2007i for a re- thefje variables may present anoma i ous i y large period change 

; ' rates (P), as a consequence of their very high evolutionary 

In the process just described, a low-mass star will typically speeds, and so at least some pre-ZAHB stars could be singled 

change its luminosity by several orders of magnitude; its tem- out in view of their high P values. As well known, high-P values 

perature may also increase very significantly. The inner struc- (e.g., |P| > 0.1 - 0.15 d/Myr) have indeed often b een observed 

ture of the star undergoes dramatic changes in the process as among field and cluster RR Lyrae stars (see, e.g., ISmith|[T995l 



1. Introduction well. And yet, in spite of the importance of this evolutionary 

phase - hereafter the pre-zero-age HB (or pre-ZAHB) phase, to 
Low-mass stars, after exhausting their central hydrogen supply, the best of our knowledge not a single such star has ever been 
develop increasingly more massive and more highly degener- positively identified. This dramatically limits the possibility of 
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for a review), whereas canonical stellar evolution theory does 
not predict such high P values for these stars, except towards the 
end of the HB phase when the star is already approaching the 
AGB (e.g., lLedll99ll) . On the other hand, small, random mix- 
ing events associated with the composition redistribution in the 
cores of HB stars have also been proposed a s possible causes for 
such h igh period changes in RR Lyrae stars (Sweigart & RenzinU 
119791) . 

The main goal of the present paper is, accordingly, to per- 
form the first systematic study of the expected properties of pre- 
ZAHB stars, with a view towards their detection using a combi- 
nation of pulsation and CMD properties. In this paper we shall, 
in particular, attempt to describe the expected pre-ZAHB stel- 
lar evolution in the case of the Galactic globular cluster M3 
(NGC 5272), which is well known for being a particularly RR 
Lyrae-rich cluster, and also for containing a number of high-P 
RR L yrae stars (see, e.g., the recent work by ICorwin & Carnevl 
1200 ll and references therein). 

In $2]we describe the method employed to compute the evo- 
lutionary tracks, along with our Monte Carlo technique to pro- 
duce synthetic pre-ZAHB distributions. In §[3] we present the em- 
pirical data that we use to compare our theoretical predictions 
with the observations. In $4]we show the results from this com- 
parison. Conclusions and final remarks are provided in $5] 

2. Methodology 

2.1. Evolutionary Tracks and Interpolation Grid 

We used the Garching Stellar Evolution Code (GARSTEC) 
to construct the evolutionary tracks required to populate the 
pre-ZAHB phase. GARSTEC is a one-dimensional, hydrostatic 
code, which does not include the effects of rotation. The program 
in its present version is capable of calculating precise solar mod- 
els as well as following the evolution of low-mass stars into the 
latest phases of their evolution, and has be en successfully us ed 
to follow a star through the He flash (e.g JSchlatfl et al.ll20Qlh . 

In order to compute the present set of evolutionary tracks, we 
adopted the mixing length theory of convection with a value of 
aMLT = 1 .75, as obtained from a c alibration of the solar model 
using the iGrevesse & Sauvall dl998l) mixture. A chemical com- 
position of X - 0.749 (where X is the hydrogen abundance) 
and Z = 0.001 was use d, as appropriate fo r M3. We used the 
OPAL equation of state (Ro gers et al.lll996|). Alexander's opac- 
ities for low temperatures ([Alexander & Ferguson 1994) , an d 
OPAL opacities for high temperatures (Iglesias & Rogers 1996). 
The observed enha ncement in the q-c apture elements among 
stars in the cluster dJohnson et al.l f2005) was taken into account 
in the comput ation of our opactit y tables . The interested reader 
is referred to ISerenelli & Weissl d2005l) and IWeiss & Schlattll 
(2000. 12007b for a much more detailed description of the nu- 
merics and input physics that go into this code. 

Our models were computed with an initial mass on the zero- 
age main sequence (ZAMS) of 0.85 M . This value was chosen 
in order to accomoda te the current estim ated age of the Universe, 
i.e. 13.7 + 0.2 Gyr dSpergel et al.ll2003l) . considering that, ac- 
cording to our computations, a 0.80 M , non-mass-losing star 
reaches the ZAHB at an age of ~ 13 Gyr. Mass l oss was applied 
in its most widely used formulation, namely the Rei mersl d 19751 
119771) formula, with a range of values for the corresponding mass 
loss "efficiency parameter" t/r that allowed us to produce from 
very red to very blue HB stars. While it is now known that the 
Reimers for mulation does not properly describe mass loss in red 
giants (e.g.. ISchroder & Cuntzl |200I ICatelanl I2008L and refer- 



ences therein), the specific choice of mass loss formalism turns 
out to be unimportant for our present purposes. 

A thin interpolation grid was created by evolving stars from 
the ZAMS to the ZAHB for a range of t]r values, from r/R = 
to ?7r = 0.6, in steps of 0.01. Following this procedure, a to- 
tal of 61 evolutionary tracks was computed. As an example of 
the results of our calculations, we show pre-ZAHB evolutionary 
tracks for different t/r values in the HRD of Figure [1] where one 
can clearly see the secondary loops associated with each succes- 
sive secondary flash (see also ©. 

According to our calculations, the total duration of the pre- 
ZAHB phase, from the RGB tip to the ZAHB, ranges from 
~ 1.12 Myr to ~ 1 .4 Myr, depending on the total mass of the star 
by the time it reaches the RGB tip. The criteria a dopted to define 
the ZAHB location are the same as proposed by Sweigart ( 1994, 
and private communication), which we now proceed to describe. 
As well known, the main He flash is produced near the center, 
but not in the exact center of the star. The He-burning region in 
the core moves inward through a series of low-amplitude sec- 
ondary flashes until it reaches the center. Associated with each 
of these secondary flashes is a temporary convective zone. This 
convection finally reaches the center of the star, but the core at 
that time is still not fully in thermal equilibrium. As a result, the 
outer edge of the convective core moves inward until thermal 
equilibrium is largely restored. A minimum in the convective 
core size is reached and the convective core starts to grow as 
expected during the HB evolution. In our calculations, we took 
the minimum in the convective core size just after convection 
reaches the center as indicating that the star has finally reached 
the ZAHB. 

We emphasize that hydrodynamical effects are ignored in 
our computations. That this is a sufficiently good ap proxima- 
tion is supp orted by the hydrodynamical calcu l ations by Deupree 
( 1996) and iDearborn. Lattanzio. & EggletorJ d2006l) . which in- 
deed show that the He flash is mainly a non-hy drodynamical 
event in the life of a low-mass star (but see also Moc ak et all 
2008, for a different viewpoint). We also note that, due to the 
fact that the He flash originally occurs off-center, a molecu- 
lar weight inversion takes place that may lead to a Rayleigh- 
Taylor instability, and thus trigger thermohaline mixing. We also 
ignore such a process in this paper, since i t has been shown 
(Kippenhahn, Ruschenp latt. & Tho mas 1980, see their §7b) that 
the timescale for such mixing to reach the center of a pre- 
ZAHB star is of order 10 s - 10 9 yr, which is much longer 
than the pre-ZAH B lifetime (~ 10 6 yr). We have repeated the 
Kippenha hn et al.l calculations in the case of our present mod- 
els, and found, on the basis of their eq. (A7), a timescale that is 
longer by about an order of magnitude than originally estimated 
by those authors. This confirms that thermohaline mixing is not 
important during pre-ZAHB evolution. 

2.2. Population Synthesis Code 

To produce synthetic pre-ZAHB distributions, a suitable mass 
distribution must be adopted. Here a Gaussian-shaped mass de- 
viate at the tip of the RGB was assumed, following the results by 
IValcarce & CatelarJ d2008l) that indicate that the mass distribu- 
tion along the M3 ZAHB does not differ markedly from a normal 
shape. Small corrections were applied to the mean and standard 
deviation values reported by those authors, due to the different 
adopted initial He content in the two studies. More specifically, 
in their calculations they used the evolutionary tracks created by 
ICatelan et all {1998) with a He content of Y = 0.23, while our 
tracks were evolved using Y = 0.25. This 0.02 difference makes 
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Fig. 1. Pre-ZAHB evolutionary tracks (solid lines), based on full evolutionary calculations (extending from the ZAMS to the ZAHB) 
using different values of the Reimers mass loss parameter rj R (indicated on the upper left side of each panel). The edges of the 
instability strip are also schematically presented (blue long-dashed lines), whereas the theoretically derived ZAHB is depicted as a 
green short-dashed line. 



our ZAHB bluer and brighter than theirs for a given ZAHB mass, 
resulting in slightly higher mass values at a given ZAHB tem- 
perature. This leads to a normal deviate with (M) = 0.647 M 
and cr = 0.022 M , as adopted in the present study. Our evolu- 
tionary track with 77R = 0.41 reaches the tip of the RGB with a 
mass value of 0.645 M , very similar to the mean value of the 
Gaussian. This track will accordingly be used in several of the 
examples to follow. 

Masses for the synthetic stars populating the pre-ZAHB 
were generated randomly following the method described in 



iBox & Mullerl d!958l) . The total pre-ZAHB lifetime for the star's 
mass is then computed by interpolating over the evolutionary 
tracks. A random number between and 1 is then generated 
for each star, with corresponding to the beginning of the pre- 
ZAHB phase (i.e., to the main He flash) and 1 to its end (i.e., 
to the ZAHB) for the track in our set whose pre-ZAHB evolu- 
tion takes longest. The star's evolutionary time (in Myr since 
the main He flash) is then computed by multiplying this random 
number by the maximum pre-ZAHB lifetime in our set of evolu- 
tionary tracks. (Underlying this procedure is the assumption that, 
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apart from random fluctuations, stars are fed into the pre-ZAHB 
phase at an essentially constant pace.) Such an evolutionary time 
is then compared against the total pre-ZAHB lifetime for the 
star's mass, as obtained previously. As a consequence, some 
stars might have an assigned age that is older than their maxi- 
mum calculated pre-ZAHB lifetime; these stars are discarded, 
since they are actually in a more evolved stage (i.e., the HB 
proper) than the pre-ZAHB phase. Once the mass of the syn- 
thetic pre-ZAHB star and its age have thus been assigned, the 
program obtains its luminosity and temperature by means of in- 
terpolations over our grid of evolutionary tracks. We are then in 
a position to check the star's variability status, and - if applica- 
ble - to compute its pulsation period. 

As can be seen from Figure Q] pre-ZAHB stars are expected 
to mostly occupy a locus similar to that of HB stars in the HRD. 
In addition, except for their innermost regions, the structure of 
these stars should be very similar to those of bona-fide HB stars 
occupying a similar position in the HRD. Therefore, any pul- 
sating pre-ZAHB stars should be very similar to bona-fide RR 
Lyrae stars, and we can accordingly use theoretical predictions 
for the instability strip edges and pulsation periods of RR Lyrae 
stars to predict the pulsation properties of pre-ZAHB stars. 

In this sense, to de fine the blue edge o f the instability strip, 
we use equation (1) in lCaputo et al.l ([1987), but applying a shift 
by -200 K to the temperature values thus derived. The width of 
the instability strip has been taken as A log T^g = 0.075; this pro- 
vides the temperature of the red edge of the instability strip for 
each star once its blue edge has been determined. These choices 
provide a good agreement with mor e recent theore tical prescrip- 
tions and the observations (see §6 in lCatelanl2004 for a detailed 
discussion). 

We include both fundamental-mode and first-overtone pre- 
ZAHB pulsators in our simulations. The computed per i ods fo r 
those stars are based on equation (4) in ICaputo et a l. (1998). 
Therefore, to compare our model prescriptions with the obser- 
vations, the observed periods of any first-overtone pre-ZAHB 
variables must be "fundamentalized," which can be achieved by 
adding 0.128 to the logarithm of the period (e.g., Catelan 2005, 
and references therein). 

Period change rates are then computed for each synthetic star 
by following in detail the variation in P with time along the in- 
terpolated evolutionary track for the star's mass. 

Let us focus on a synthetic star with mass M = 0.648 M , 
and age since the He flash randomly assigned as 1.119 Myr. By 
interpolation on our grid, we find for this synthetic star a position 
in the HRD characterized by log(L/L ) = 1.711, log(r e ff) = 
3.827. This places the star inside the instability strip, pulsating 
with a period P = 0.546 days - the latter decreasing at a rate 
P = -0.389 d/Myr. 

The result of the interpolation procedure is shown in 
Figure [2] According to the star's mass value, it should be po- 
sitioned between our evolutionary tracks for 77R = 0.40 and 0.41 . 
The figure shows the time evolution of the physical parameters 
along these tracks (blue for 77R = 0.40 and red for 77R = 0.41). 
The time range has been selected to emphasize the time when 
the tracks cross the instability strip (i.e., 3.75 < log T e s < 3.85, 
approximately). As can clearly be seen, the star's position (in- 
dicated as a black dot) has been correctly interpolated between 
the two tracks. From the third panel, it can also be appreciated 
that changes in the trend of increasing or decreasing period val- 
ues are indeed accompanied by zero values for P (fourth panel), 
and also by changes in the slope of the temperature curve (sec- 
ond panel). These changes correspond to the secondary flashes, 
which produce the previously noted wide loops along the evolu- 
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Fig. 2. Time evolution for the 77R = 0.41 (red, solid) and 77R = 
0.40 (blue, dashed) tracks. The black dot represents the position 
of one of the stars, as generated by our interpolation code, whose 
mass places it between these two evolutionary tracks. Upper 
panel: luminosity evolution. Second panel: temperature evolu- 
tion. Third panel: period evolution. Lower panel: period change 
rate evolution. See text for details. 



tionary tracks in the HRD. It is also interesting to note that, at 
least in principle, very high (negative) period change rates (e.g., 
of order -15 days/Myr) are possible, although not very likely 
(see Fig. |5]l - and indeed, no stars with such extreme P values 
have been detected in M3. We will come back to this point later 
in this paper. 



3. Empirical Data 

To effectively compare our theoretical predictions with the ob- 
servations, we need a database containing both variable and non- 
variable stars. Mag nitudes and colors for the variable stars in M3 
were obtained from lCacciari et all (120051) . who provided quanti- 
ties for the "equivalent static star" for a total of 133 variable stars 
in the cluste r. For the non- variable s tars, we used the photom etry 
provided bv lBuonanno et al.l d 19941) and lFerraro et al.l dl997l) . 

Our evolutionary tracks were converted from the theoreti- 
cal [log(L/L Q ), log(r e ff)] to the observational [My, (B-V)] plane 
using th e color transformations and bolometric corrections pro- 
vided bv lVandenBerg & Cleml (120031) for a metallicity [Fe/H] = 
-1.5 and abundance of the alpha elements [a/Fe] = +0.3. 
In terpo l ations on their tables was carried out using the algorithm 
bv lHilll d 1982b . We converted the thus derived magnitudes to ap- 
parent magnitudes using a distance modulus of 15.14 mag in V, 
based on a comparison between our evolutionary tracks and the 
observation al data. This i s only slightly different from the value 
provided by lHarrisI d 1996b . namely 15.12 mag. Magnitudes and 
colors were corre cted for extin ction based on a reddening value 
E(B-V) = 0.010 dHarrisll 19961) . and adopting a standard extinc- 
tion law with A v /E(B-V) = 3.1. 
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Fig. 3. Results for one simulation with 457 pre-ZAHB stars. Red 
triangles indicate the synthetic pre-ZAHB stars. Upper Panel: 
The solid line depicts the evolutionary track for a value of the 
Reimers mass loss parameter of t/r = 0.41, whereas the green 
dotted line shows the position of the ZAHB for different mass 
values. Lower Panel: As in the upper panel, but zooming in 
around the secondary loops of the evolutionary track, and in- 
cluding the instability strip boundaries (blue dashed lines). 



4. Results 

M3 has ~ 530 HB stars w hose positions on the CMD have been 
measured dCatelanll2004l and references therein). As shown in 
Figure[3] the pre-ZAHB variable stars are expected to be primar- 
ily located above the ZAHB, and thus at least partially superim- 
posed on the locus occupied by the bona-fide HB stars (co mpare, 
e.g., with the CMD presented bv lCorwin & Carnevll200ll) . 

If we consider lifetimes, pre-ZAHB evolutions lasts, for the 
range of masses we are interested in, ^1.3 Myr. We evolved 
our 77r = 0.41 star up to the end of the HB phase, where He 
is exhausted in the core, and the evolution took 77 Myr. 
Consequently, one should expect to find one pre-ZAHB star for 
every 60 or so bona-fide HB stars. Therefore, among the 530 HB 
stars in M3, approximately 9 should be in the pre-ZAHB phase - 
two of which being pulsating (RR Lyrae-like) stars. Do the lat- 
ter's expected pulsation properties differ in any noteworthy way 
from those characterizing bona-fide RR Lyrae stars? 

4.1. Expected Pulsation Properties of Pre-ZAHB Stars 

To answer this question, we ran pre-ZAHB Monte Carlo simula- 
tions for an M3-like input mass distribution. We used a total of 
600 stars, of which ~ 450 are still in the pre-ZAHB phase, the re- 
maining having already evolved past the ZAHB. (Since the age 
of the star is generated randomly, some of them do fall in the 
pre-ZAHB phase for a given mass value, while others have al- 
ready left this evolutionary stage.) The latter are discarded, since 
we are only interested in the evolution prior to the ZAHB. The 
output from one of these simulations is shown in the HRD of 



Fig. 4. Distribution of fundamentalized periods from our simu- 
lations, for the 1905 pre-ZAHB variable stars obtained from our 
simulation for a total of 9427 pre-ZAHB stars 



Figure[3] where the position of the instability strip is also shown, 
along with the pre-ZAHB track for 77R = 0.41, as well as the 
ZAHB locus obtained from our tracks for different mass values. 

As can be seen from this plot, there is a considerable number 
of pre-ZAHB stars lying within the boundaries of the instability 
strip, and which can accordingly be tagged as variable stars. In 
this simulation approximately 100 variable stars were produced, 
corresponding to 22% of the total number of pre-ZAHB stars in 
the simulation. Therefore, in the case of a globular cluster like 
M3, we expect that approximately one in every four or five pre- 
ZAHB stars will be disguised as bona-fide RR Lyrae stars. 

In Figures |4] and [5] we present histograms that show the ex- 
pected period and period change rate distributions for an M3- 
like input mass distribution. To minimize the effects of statisti- 
cal fluctuations, we increased the number of stars in the Monte 
Carlo simulations up to a total of 10000 stars, 9427 of which 
turned out to be on the pre-ZAHB phase. The latter are the stars 
that are plotted in these figures. 

It is worth noting that, according to the derived period his- 
togram, the expected periods of pre-ZAHB pulsators are rela- 
tivel y long, especially i n view of the fact that M3 is the prototype 
of an Oosterhofn (1 1 9391) type I cluster, where the average periods 
of th e fundamental-type RR Lyrae are of order 0.55 days (see, 
e.g.. lClement et al]|200ll) . Indeed, the peak of the fundamental- 
ized period distribution is found at 0.55 d in M3, and there are 
only a few RR Lyrae stars in the cluster wit h periods long er than 
0.65 d (see the corresponding histogram in Catelan 2004). 

Period changes are related to changes in the structure of a 
star through the period-mean density relation. More specifically, 
the periods of variable stars whose temperatures are decreasing 
in time, or whose luminosities are increasing, are expected to 
also increase. These concepts successfully explain the average 
behavior observed in globular clusters, where objects with bluer 
HB types - where a higher number of redward-evolving RR 
Lyrae stars is expected - tend to have higher average P values 
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Table 1. Average Colors, Magnitudes, and Periods as a Function of P 



Range (d/Myr) 


<V) 


<T(V) 


(B-V) 


o~b-v 


<P)(day) 


o>(day) 


P < -0.8 


15.56 


0.054 


0.343 


0.062 


0.70 


0.14 


-0.8 < P < 


15.66 


0.015 


0.320 


0.061 


0.59 


0.11 


P < 


15.61 


0.061 


0.331 


0.063 


0.64 


0.14 


P>0 


15.63 


0.053 


0.338 


0.061 


0.66 


0.13 



3.5 




Fig. 5. As in Figure HI but for the period change rate distribu- 
tion. The middle panel presents the same data as the upper panel, 
but zooming in around the peak of the distribution, whereas the 
lower panel shows the wing of the distribution for the more ex- 
treme (negative) P values. 



than t hose with redder HB morphology (see, e.g., §9 in ICatelanl 
2005). As can be seen from Figure[5] most pre-ZAHB stars have 
period change rates contained in the [-1,0] d/Myr range, mean- 
ing that evolution occurs mainly towards higher temperatures 
and/or decreasing luminosity - as expected for a pre-ZAHB star 
(see also Fig. U}. 

To better understand this result, the time distribution along 
the evolutionary track corresponding closely to the most likely 
mass among M3 pre-ZAHB stars is presented in Figure [6] In 
this figure, a red dot was placed for every 0.01 Myr of evolution, 
starting from the He flash and up to the ZAHB. As we have seen, 
the period changes along the blue loops can reach values up to 

15 d/Myr, which is due to the very fast evolution of stars in 

the secondary loops. However, this is a relatively unlikely situa- 
tion (see the lower panel in Fig. [5), since pre-ZAHB stars spend 
most of their time in the blueward evolution that takes place 
after the last secondary flash, immediately before reaching the 
ZAHB. This implies that most (i.e., 76%) pre-ZAHB variables 
will actually have negative P values - as confirmed by Figure [5] 
which in addition reveals that the most likely P value for a pre- 
ZAHB star is « -0.3 d/Myr. As a consequence, according to 
our simulations, only 24% of the pre-ZAHB stars, for an M3- 
like HB morphology, are expected to present positive P values. 
On the other hand, 38% of the pre-ZAHB pulsators present P 



Fig. 6. Upper Panel: Pre-ZAHB evolution on the HRD for the 
star with 77R = 0.41. The evolutionary track is shown as a solid 
line, and red dots are plotted for every 0.01 Myr of evolution 
since the main flash. Lower Panel: As above, but zooming in 
around the secondary loops, and including the instability strip 
edges (dashed lines). 



values more negative than -0.8 d/Myr. We again emphasize that 
such extreme P values are completely unexpected on the basis of 
canonical HB evolution, according to which period change rates 
should not be higher (in absolute value) than abo ut 0. 1 d/Myr in 
a cluster such as M3 (see, e.g., Fig. 2 in Lee 1991). 

A closer look at our results is provided in Table [TJ where we 
show the average V-band magnitudes and B-V colors, along with 
the mean fundamentalized periods (and corresponding standard 
deviations), for different groups of stars with different P values 
in our M3 simulations. The most noteworthy trend present in 
this table is for the stars with P < -0.8 d/Myr to be slightly 
brighter, and have correspondingly longer periods, than those 
with smaller P values. However, these trends are far from consti- 
tuting firm tendencies, given the fairly large standard deviations 
that accompany the studied quantities. 

In Figure [7] we present a plot s howing the variable stars in 
M3 (data from ICacciari et alj 120051) . together with the ZAHB 
and the pre-ZAHB evolutionary track for t]r = 0.41. Stars with 
measured period change rates higher than +0.15 d/Myr are plot- 
ted as open circles, whereas those with P < -0.15 d/Myr are 
shown as filled circles. In Figure[8]we present the same variable 
stars but now including three evolutionary tracks: for 77R = 0.42, 
77R = 0.41, and 77R = 0.40. These tracks were selected because 
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Fig. 7. Upper Panel: CMD for M3, including only the RR Lyrae 
variable stars. Filled red circles: stars with P < -0.15 d/Myr; 
open blue circles: stars with P > +0.15 d/Myr; triangles: RR 
Lyrae stars with smaller (or undetermined) P values. The evolu- 
tionary track (solid line) for an t]r = 0.41, as well as the ZAHB 
(green dashed line), are overplotted on the observed data. Lower 
Panel: Same as in the previous panel, but zooming in around the 
instability strip. 



their associated mass values at the tip of the RGB (0.639 M , 
0.645 M Q , and 0.650 M Q , respectively) make them representative 
of the majority of pre-ZAHB stars in M3, whose HB mass dis- 
tribution is characterized by (M) = 0.647 M and <x = 0.022 M 
(see 32.2l for details). As can be seen from both these figures, the 
positions of these variables in the CMD are broadly consistent 
with at least some of them being pre-ZAHB stars - and indeed, 
according to our simulations, at least ~ 2 of them very likely are. 

4.2. Variable vs. Non-Variable Pre-ZAHB Stars 

As we have seen, our simulations indicate that of order 20% of 
the pre-ZAHB stars are expected to be RR Lyrae-like variables. 
What about the remaining 80% of these stars? 

To answer this question, we divide the pre-ZAHB stars into 
four different groups, depending on their photometric and pulsa- 
tion characteristics, as follows: RR Lyrae-like, red HB-like, blue 
HB-like, and AGB/RGB like. The relative proportions of pre- 
ZAHB stars in each of these groups, according to our simulation 
with 9427 pre-ZAHB stars, are provided in Table [2] 

According to Table |2l most pre-ZAHB stars are expected 
to fall in the red HB region, with a proportion of almost three 
red HB-like pre-ZAHB stars for every RR Lyrae-like pre-ZAHB 
star. Comparable fractions of blue HB-like and RR Lyrae-like 
pre-ZAHB stars are also implied. But a minor fraction (1.4%) is 
expected to be sufficiently red and bright (i.e., brighter than the 
expected terminal-age HB, or core helium exhaustion, locus) as 
to fall into the AGB/RGB-like category. This is a consequence of 
the extremely fast evolutionary rates immediately after the pri- 
mary He flash (see Fig. [6]). In the case of M3, therefore, taking 
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Fig. 8. Distribution of M3 RR Lyrae stars on the CMD (symbols 
are the same as in Fig. with three evolutionary tracks super- 
imposed: dashed magenta line: track for ^r = 0.42; solid line: 
track for r]Yt, = 0.41; dash-dotted blue line: track for ^r = 0.40. 
Masses at the tip of the RGB are 0.639 M , 0.645 M , and 
0.650 M Q , respectively. See text for mode details. 

Table 2. Relative Proportions of Pre-ZAHB Stars in M3 
Simulations with 9427 Pre-ZAHB Stars 



Classification 


Number 


Number fraction 


Red HB-Like 


5407 


57.4% 


RR Lyrae-Like 


1905 


20.2% 


Blue HB-Like 


1979 


21.0% 


AGB/RGB -Like 


136 


1.4% 



into account a total population of 530 HB stars, it appears most 
likely that there are ~ 2 RR Lyrae-like, ~ 2 blue HB-like, ~ 5 red 
HB-like, and an insignificant number - most likely 1 at most - 
AGB/RGB-like stars. 

5. Conclusions 

In this paper, we have performed a detailed theoretical study of 
the expected evolutionary and pulsation properties of pre-ZAHB 
stars for a metallicity appropriate to the case of the Galactic 
globular cluster M3. Our results indicate that there should be 
around one pre-ZAHB star for every 60 or so bona-fide HB stars 
in the cluster. Based on extensive Monte Carlo simulations we 
find that, among the pre-ZAHB stars, of order 20% should ex- 
hibit RR Lyrae-like pulsations - to be compared with 21% and 
57% blue HB- and red HB-like pre-ZAHB stars, respectively, 
and but a very small (of order 1 - 2%) population of AGB/RGB- 
like pre-ZAHB stars. 

The pre-ZAHB, RR Lyrae-like pulsators are characterized by 
relatively long periods (typically longer than the average for the 
cluster), and (especially) by period change rates falling mostly in 
the [-1,0] d/Myr range, with -0.3 d/Myr being the most likely 
value. However, more extreme P values are also possible; in fact, 
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38% of the variables in our simulation have P < -0.8 d/Myr. 
Relatively few (i.e., around 24%) of the pre-ZAHB pulsators are 
found with positive period change rates. 

Our pre-ZAHB Monte Carlo simulations show, in addition, 
that for an M3-like HB mass distribution the pre-ZAHB CMD 
positions should be basically indistinguishable from those of 
bona-fide HB stars in the cluster. We predict that, in the specific 
case of M3, ~ 9 pre-ZAHB stars should be present (assuming a 
total population of 530 HB stars), ~ 2 of which most likely dis- 
guised as RR Lyrae stars with relatively long (and negative) pe- 
riod change rates. Since t he number of stars in M3 with high P is 
significantly larger (e.g.. lCorwin &Carnevl l2001). we conclude 
that pre-ZAHB evolution may explain the high-P phenomenon 
for some (but not all) of the RR Lyrae stars with high period 
change rates in this cluster. 

In the future, we will extend our calculations to other clusters 
(and dwarf galaxies) with different metallicities and HB mor- 
phologies. These calculations should also prove important for 
the identification of candidate pre-ZAHB pulsators among field 
stars in the Galaxy and neighboring galaxies. With the era of 
the survey telescopes quickly approaching, the number of such 
stars available for study should dramatically increase, thus mak- 
ing empirical studies of pre-ZAHB evolution increasingly more 
feasible - at least in the long term. 
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